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Abstract 

Long-time behavior of spatial power-law correlation in one-dimensional 
self-gravitating system is numerically investigated. The power-law 
structure persists even after the system is virialized. The structure 
gradually disappears with energy exchange among particles. Lifetime 
of the power-law structure is estimated to be proportional to the sys- 
tem size. 
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1 Introduction 

Spatial structure in many-body systems is an important subject in many 
branches of physics, such as molecular systems, gravitational systems, and 
so on. For various types of such structures, their origins are interesting 
subject of study. Some of them may be realized as equilibrium distribution, 
and others may be already set at initial conditions. 

Recently we have discovered that spatial structure with power-law cor- 
relation spontaneously emerges from uniformly random initial conditions in 
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one-dimensional self-gravitating system What is important in this novel 
phenomenon is that the spatial structure is not given at the initial condition, 
but dynamically created from a state without spatial correlation. Succeeding 
research clarified that the structure is created first in small spatial scale then 
grows up to large scale through hierarchical clustering ||. 

The question to be discussed now is the properties of the state which 
has power-law correlation. In particular we have to clarify whether the state 
persists for long time, or just a short transient and disappears quickly. If the 
state lasts for long time, we have good chance to observe it and it has much 
importance. On the other hand, if the state disappears quickly, it may not 
be observable easily and it is less important to the whole dynamics of the 
system. 

Another important question is on the relaxation process ( See [§-0] and 
references therein.). Since the system is an isolated many-body system with- 
out dissipation, one would expect that the system will eventually relax to 
thermal equilibrium for sufficiently long time scale. Hence relation between 
decay (if it does) of power-law structure and relaxation processes of thermal- 
ization is also important. 

The purpose of this paper is to investigate the long-time behavior of the 
power-law structure; its lifetime and mechanism of relaxation. The model 
is introduced in the next section. In section [3] we introduce several quanti- 
ties which represent relaxation and investigate their dynamics to clarify the 
mechanism by which the power-law structure disappears. The final section 
is devoted to summary and discussions. 
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2 Model 



The model we use in this paper is the self-gravitating one-dimensional sys- 
tem 0-0. It represents a system of particles on one-dimension interacting 
through Newtonian gravity. It can also be interpreted as a system consists 
of equivalent mass-sheets placed parallel in 3-dimensional space, hence it is 
also called as the sheet model. 

The Hamiltonian of the model is 

N p 2 

H = — — + 2irGm 2 \xj — Xj \ , — oo < Xi < oo , (1) 

i=l i>j 

where X{ and pi represent coordinate and momentum of z'th particle, N is the 
number of particles, m stands for mass of each particle, and G is a positive 
constant. Force between two particles is attractive and long-ranged. 

In this paper we set m = 1/N and AttG = 1 . Time is measured in unit 
t c = l/ \JAtiGM / L , where M = mN is total mass, L is the spatial length on 
which particles are distributed initially. 

3 Numerical results 

Now we show temporal evolution of the power-law correlation. Since the 
formation process is already studied in our previous work ||, here we focus 
on the long-time behavior. 

In what follows we use a typical example whose initial condition and 
parameters are defined as follows: 

system size : iV = 2048 , 
initial condition : Wi Vi = 0, (2) 
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Xi = uniformly random in [0, 1). 
3.1 Power-law correlation 

In our previous works we showed that spatial structure where two-point cor- 
relation function obey power-law 

£(r) ex r- a (3) 

is formed in the model ([!]) from uniformly-random initial conditions II], 0- 
Here the the two-point correlation function £ (r) is defined as 

dP = ndV{\ + £(r)) (4) 

where dP stands for probability to find another particle in volume dV at 
distance r from a particle, n is the average number density. 

Using this correlation function we found that the model (P with initial 
conditions with zero velocity dispersion evolves into a state where £(r) obeys 
power-law (|3]). In Figure [l| we show a typical example of such power-law 
correlation. 




Figure 1: An example of power-law behavior of correlation function (f|) for 
parameters t = 250. £(r) oc r -0 ' 11 . 
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3.2 V r : virial ratio 

The main question we would like to discuss in this paper is whether the 
power-law structure lasts for long time or not. For that purpose we need to 
define a time scale for reference. 

As for the one-dimensional self-gravitating system ([[]), it is known that 
relaxation proceeds through three stages; first virialization process occurs, 
then "microscopic relaxation", then "macroscopic relaxation" f4j,|J. (Each 
process is explained later.) We will compare the decay of power- law structure 
and each relaxation process. 

First we adopt a time scale called virial time t vr defined from virial ratio 
V r (t). Virial ratio is defined as 



where Eki n {t) and E pot (t) are kinetic and potential energy of the system at 
time t, respectively. 

Since the total energy is conserved and finite, it follows from virial theo- 
rem that 



where E kin and E pot are long time average of E kin {t) and E pot (t), respectively. 
Thus convergence of V r (t) to 1 is one good measure of relaxation of the 
system, and the time scale t vr at which V r (t) converges to 1 is a good reference 
timescale. We call t vr virial time. 

Figure shows temporal evolution of the virial ratio It is seen that 
the virial ratio V r (t) converges to 1 at about t ~ 150. At this time it can be 
said that the system has relaxed to a state in the sense of energy partition. 



V r (t) = 2E kin {t) / E pot {t) 



(5) 




(6) 
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Figure 2: Time evolution of virial ratio for parameters (0). 

The important point to note is that the power-law correlation shown in Fig.|l] 
exists after the virial ratio is converged to 1. (Remember that the figure is 
taken at t = 250. Also it is confirmed in other examples.) That is, power- 
law correlation lasts longer than virial time. In this sense the power-law 
correlation persists long. 

3.3 a(t) : exponent of power-law 

The next problem is the length of time it lasts. For that purpose we investi- 
gate the long-time behavior of the exponent a of power- law 

Fig. |3] represents temporal evolution of the exponent a for time scale 
much longer than the one shown in Fig. |2| for 10 independent orbits. Initial 
conditions are the same as (0) with different realization of random numbers. 
From this figure we see that the power-law structure is not stationary but 
transient and gradually fading in this time scale. Decrease of the exponent 
a(t) can be considered as relaxation process. 
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Figure 3: Long-time behavior of the exponent a(t) of correlation (§) for 
< t < 10000 for 10 initial conditions. 
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Figure 4: Temporal evolution of A(t) for < t < 10000 . 
3.4 A(t) : energy exchange 

Next we investigate the process of relaxation of power-law correlation. For 
that purpose we measure A(t) defined below, by which we can see the degree 
of energy exchange between particles f|,||. 

A(t) is defined as follows. First we define £j(i), 1-particle energy per unit 
mass as 

1 N 
e i (t) = -v 2 i {t)+2T:GmY J \x J {t)-x i {t)\ . (7) 

1 3=1 

If the system is ergodic the long time average of £i(t) converges to a 
unique value 

lim i f 7 £l (t)dt = 5E total /3 = e (8) 

T-*oo 1 Jo 
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for all i. 



Using this e , A(t) is denned as 



l l N n I* 

A W = -\ T?E T / e t {t')dt'-e, 
eq \ N fr[Vt Jo 



(9) 



A{t) represents deviation from equipartition of energy among particles. At 
each time step Ei(t) varies from particle to particle. A particle which initially 
has high energy may give the energy to particles which initially have low 
energy, and vice versa. Hence, in the course of time evolution, if the value 



\t J £i ^'^^' ~ £ °j decays t° zero f° r m °st of the particles, we can say that 
energy is well exchanged among particles. On the other hand, A(t) remains 
constant if energy is not exchanged between particles and every Si{t) is kept 
constant. Thus decrease of A{t) indicates energy exchange among particles. 

Fig. H shows temporal evolution of A(t) for the parameters Initial 
increase of A(t) represents that, instead of becoming equipartitioned, differ- 
ences of one-particle energy among particles are temporarily enhanced. Let 
us turn our attention to behavior of A(t) in later time. We see that A(t) 
begins to decrease from t ~ 1500, which is of the same order of the time 
when the power-law fades away as seen in Fig. |3|. Hence we can say that the 
power-law structure relaxes as energy is exchanged among particles. 

3.5 v{e) : cumulative energy distribution 

Now we investigate further the process of relaxation of power-law structure. 
We have seen that the power-law structure begins to disappear as energy 
exchange among particle begins, as indicated by the decrease of A(t). If, 
through this relaxation, energy distribution is kept constant, this relaxation 
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phase corresponds to " microscopic relaxation" discussed by Tsuchiya et al. 
J|,[§]. If, on the other hand, energy distribution is also relaxed to that of 
thermal equilibrium, then the phase corresponds to "macroscopic relaxation" . 
We measure energy distribution by v(e) defined as follows: 

u(e) = — (number of particles whose energy is less than e) , (10) 

where N is the total number of particles. Since u(e) is a cumulative distri- 
bution, v(e) increases monotonically with respect to e. If the graph of v(e) 
does not change with respect to time through relaxation, then the relaxation 
is microscopic relaxation. 



cumulative 1 -particle energy distribution 




Figure 5: v{e) for t = 250, 300, • • • , 2000 



cumulative 1 -particle energy distribution 




1 -particle energy 



Figure 6: v{e) for t = 250, 300, • • • , 2000, 8000, 8050, • • • , 9750 . 

Figs. [5] and |6] show u(e) for two temporal intervals. Two well-known quasi- 
stationary state, water-bag distribution and isothermal distribution J3],[5|] are 
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also shown as references. 

Fig. [5] represents u(e) for 250 < t < 2000, when the system has power-law 
correlation and the system is virialized (see Fig. |2|). Please note that u(e) 
does not change for this time. 

Fig. H represents v(e) for 8000 < t < 9750 overlaid on the previous 
figure. In this temporal interval power-law correlation is lost (see Fig. It 
is clearly seen that the energy distribution remains the same regardless of 
the existence of power-law correlation. Thus we can say that the relaxation 
process of power-law structure can be explained by the microscopic relaxation 
caused by energy exchange among particles while keeping the global energy 
distribution fixed. 

4 Summary and discussions 

In this paper we numerically investigated the long-time behavior of power- 
law correlation structure |l|, |2|] in one-dimensional self-gravitating system. 
The power-law correlation exists even after the system is virialized. In this 
sense the structure persists for long time. The lifetime of the structure is, 
however, not infinite and it gradually fades away. By comparing the timescale 
of the lifetime with other time scales of relaxation we see that the decay of 
the power-law structure is not caused by virialization process but by energy 
exchange among particles. The one-particle energy distribution v(e) is in- 
variant through the relaxation of the structure, hence it is the microscopic 
relaxation H,|J. The system is not stationary if we see it by the exponent a 
of correlation, but it is at the same time quasi- stationary if we see it by the 
energy distribution v(e) after the system is virialized. 
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The fact that power-law correlation structure lasts long implies that the 
structure is important in the whole dynamical evolution of the system and 
we may be able to observe such structure in other systems even if the system 
is virialized. 

Relation between the decay of power-law correlation and microscopic re- 
laxation has an important consequence. It is known that the microscopic 
relaxation time depends linearly on the system size t micro oc N ||. Hence 
it is likely that the average lifetime of the power-law correlation has similar 
dependence on N, which is quite long for large N. 

Let us now summarize schematically the entire life of the spatial power- 
law correlation structure in the model (0) clarified through a series of our 
works. 

1. From initial condition with zero velocity dispersion |]J, spatial struc- 
ture with power-law correlation emerges in small spatial scale 0. Typ- 
ical example of such initial condition is Vi v $ = , Xi = random. 

2. The power-law correlation structure grows up to large scale through 
hierarchical clustering ||. 

3. The system is virialized. Power-law correlation structure persists. 

4. Microscopic relaxation takes place and power-law correlation structure 
disappears. Energy distribution remains invariant. 



Related models with a uniform background or friction [ pl2] , |T3|| are inves- 
tigated from cosmological interest. In those models similar fractal structure 
has been observed and is more persistent. Comparison with our results may 
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be helpful to understand mechanism by which the structure is retained in 
their models. 

It should also be noted that the structure formation we have studied have 
origin in dynamics. It is interesting that, instead of monotonously relaxing to 
thermal equilibrium the system, initially does not have spatial correlation, 
creates spatial correlation of power-law type through dynamics. It seems 
difficult to explain it in terms of thermodynamics or equilibrium statistical 
physics, since the structure appears in states which are not thermally relaxed. 
Recently dynamical formation of spatial structure in long-range interacting 
systems is actively investigated. Relation between our results and other 



examples of structure formation with dynamical origin |[TQ|, |ll],[i^,|I5|j will be 
interesting. 

In this paper we have focused on typical behavior of the system. Several 
problems related to details of the behavior remain unsolved. Among them 
there are, for example, initial condition dependence, system size dependence, 
and so on. These are beyond the scope of this short paper and of important 
subjects of future study. 

In this paper we have studied the behavior of the system until power-law 
structure disappears through energy exchange and energy distribution. It 
will be interesting to investigate whether there are other relaxation processes 
effectively working ||T6|,|T7|1, and the behavior of the system in the course of 
approach to equilibrium []|,|,[T8|,[19 . 



The spontaneous emergence of power-law correlation we have studied so 
far is a novel phenomenon with rich implication. We hope our research serves 
as a foundation for future study, both fundamental and applied. 
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